Ergodic Considerations in the gravitational potential of the Milky Way 
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ABSTRACT 

A method is proposed for constraining the Galactic gravitational potential from high precision observations 
of the phase space coordinates of a system of relaxed tracers. The method relies on an "ergodic" assumption 
that the observations are representative of the state of the system at any other time. The observed coordinates 
serve as initial conditions for moving the tracers forward in time in an assumed model for the gravitational 
field. The validity of the model is assessed by the statistical equivalence between the observations and the 
distribution of tracers at randomly selected times. The applicability of this ergodic method is not restricted 
by any assumption on the form or symmetry of the potential. However, it requires high recision observations 
as those that will be obtained from missions like SIM and GAIA. 

Subject headings: Astrometry — Galaxy: halo — Galaxy: kinematics and dynamics 



1. Introduction 

Measurements of velocities of cosmological objects are 
a classic probe of the mass distribution on all scales. It 
was the motions of individual galaxies in galaxy clusters 
which first showed that luminous mater contribu ted only 
a sm all fraction to the total mass in clusters ( Zwickvl 
Il937[ ). implying the existence of dark matter. The com- 
bined mass of the Milky Way (or the Galaxy) and M31 
is constrained from the observed relative motion between 
the two galaxies and the requirement that the initial dis- 
tance between thei r respective centers of mass vanishes 
near the Big Bang (jKahn fc Woltjeij llQSQY Line-of-sight 
velocities of other galaxies in the Local Group of galax- 
ies also are used to estimat e its mass by the condition of 
vanishing initial distances (jPeebled'lOSQ"). On scales 10s 
of Mpcs, peculiar motions (deviations of Hubble flow) of 
gala xies constrain the global mass density in the Universe 
re.g. lNussedl2nn8h . 

Determining the mass distribution in the Milky Way 
is particularly important. There is ample information 
on the baryonic content of the Galaxy which could be 
modeled in detail only if the dark matter distribution 
is known. The rotation curve of the Galaxy is limited 
to distances smaller than 20 kpc and does not provide 
any information on deviations from spherical symmetry 
of the halo. The mass distribution at larger distances, 
motions of Galactic satelli te galaxies, globular clusters 
and stars are invoked (e.g. ISakamoto et al ][200i. Con- 



straints on the Galactic mass from these tracers are de- 
rived from the condition that the observed speeds of 
Galactic objects do not exceed the escape velocity. This 
approach yields only a lower limit and is mostly sensi- 
tive to the highest velocity objects. The vast majority 
of the sample objects play no role in deriving the mass 
limit. Alternatively, one could adopt a Bayesian like- 
lihood formalism in which the phase space distribution 
function is assumed to follow certain form which could 
be matched with th e observations to probe the Galactic 
potential field (e.g. 'Little fc T remaing Il987t iKochane^ 
[1996; Wilkinson & Evans 1999) . 

Proper motions, resulting from velocities perpendic- 
ular to the line-of-sight , are curre ntly measured only 
for nearby tracers (e.g. ISakamoto~e t al. 2003). There- 
fore, most current mass estimates rely on the measured 
line-of-sight motions of individual components. During 
the next few years, accurate proper motions for a large 
sample of Galactic tracers are expected to be measured 
by the space missions Glo bal Astrometry Interferometer 
for Astrophysics (GAIA) (jLindegren fc PerrvmanI 19961) 
and Space Interferometry Mission fSIM) (^ Unwin et al.l 
|2008[ ). Even accurate phase space information require 
additional assumptions in order to constrain the Gravi- 
tational potential 0. We present here a general method 



^ The gravitational force field is equal to the acceleration rather than 
velocities. Measuring the acceleration of a tracer with orbital pe- 
riod over an observing time t^^^ requires astrometry with angu- 
lar resolution a factor t^j^^ /i ^ higher than the precision needed for 
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which rehes on high precision measurements of positions 
and velocities of tracers. The method assumes that the 
tracers have reached an equihbrium state in the Galactic 

gravitational field. 

2. The method 

The expected precision of future data of proper mo- 
tions, radial velocities and distances will allow an ac- 
curate determination of the orbits of tracers in a given 
Galactic gravitational field. This motivates the "ergodic" 
method which constrains the gravitational potential by 
the hypothesis that the state of a dynamically relaxed 
system of tracers at any time is statistically representa- 
tive of that system at any other time. The method is 
described as follows. 

1. assume a model for the gravitational potential field 
^'(r, a, /3 . . .) where a, /? etc are free parameters. 

2. for a given choice of free parameters, use the ob- 
servations as initial conditions at time i = to 
advance the particles (tracers), in the gravitational 
field ^, for sufficiently long time. 

3. select snapshots of the particle distribution at times 

t>t„ 

4. compare between statistical measures of the parti- 
cle distributions in the observed data (t = t^) and 
in the snapshots {t > t^). If necessary, repeat (ii)- 
(iv) with a different choice of free parameters until 
a reasonable agreement is reached. 

We work here with four statistical measures. The first is 

the N statistic computed as follows. Consider all snap- 
shots (at t > tg) corresponding to an assumed value of a. 
For a snapshot s, we compute A^s defined as the number 
of particles with Galactic distances that are larger than 
the respective observed distances. We define the statistic 
N as N =< Ng >, the average of Ng over all snapshots. 
We further define a =< (Ng — iV)^ >^/^, as the r.m.s. 
scatter in Ng. In the limit of a large number of tracers, n, 
the limits = n/2 and N = n/2 are approached for the 
correct model. The second is the distribution function of 
the Galactic distances, r, of tracers, i.e. the density of 
tracers as a function of distance. The third measure is 
distribution function of 



velocities. Since t^j^^ ~ a few years while ~ a few Gyrs, the task 
is out of reach in the near future. 



computed from the Galactic distances, r, of particles in 
each snapshot. The pericenter, r^.^, and apocenter, r^^^, 
of each particle is computed from the numerically inte- 
grated orbits. The distributions of r and ^ computed 
from the observations given at tg should be statistically 
equivalent to the respective distributions in any snapshot 
a.t t > tg if the system is evolved to t > tg using a gravi- 
tational potential which is a reasonable approximation to 
the true potential field. Significant differences between 
the initial and later distributions should appear if the 
integration at t > tg is performed with sufficient devi- 
ations from the original gravitational potential. To see 
this consider the simple case of particles executing circu- 
lar motions in the gravitational field of a point mass, Mq. 
Now let us use a snapshot of the motion of this system at 
time t = tg as initial conditions to move the particles for- 
ward in time in the field of the point mass M ^ Mq. The 
pericenters, r^.^ , of orbits obtained with M will equal the 
initial distances so that ^ = for all particles at time tg . 
However, the distances at a randomly selected later time 
t ^ tg span the range r^.^ to r^^^, yielding ^ from to 
1 (the ^ distribution peaks at the end points where the 
radial velocities are zero). Therefore, M <C Mq could 
be rejected by the gross mismatch between the ^ distri- 
butions at the initial time tg and at i ^ tg . Similar 
arguments apply for the distribution function of r. 

The fourth statistic is based on the virial relation 

between the total kinetic energy (per unit mass) T = 
^r^/2 and a term W related to the gravitational force 
field. In order to derive the explicit relation, we per- 
form a scalar multiplication of the radius vecor, r of 
each particles with the corresponding equation of mo- 
tion dv/dt = — V^*. Summation over all particles and 
time averaging gives 

2T = W, where VF=- ^ r-V*, (2) 

tracers 

and the over-lines denote time averaging. Therefore, as 
our fourth statistic we consider the ratio \W\/2T com- 
puted directly from the observations for an assumed form 
of "if. Since the virial relation applies to time average 
quantity, this statistic, computed at a single time, is ex- 
pected to deviate from unity for a finite number of parti- 
cles. Therefore, although the statistic could be computed 
without evolving the particles in time, snapshots at later 
times will be used to compute the scatter in the statistic. 
This ratio differs from the other statistics in that it di- 
rectly involves the observed velocities. It is only sensitive 
to the force field at the observed positions of the tracers, 
i 



2 



3. Tests 

Thorough tests of the performance of the method 
should take into account the details of future observa- 
tions. Although this is necessary in order to provide 
reliable error estimates, the task seems futile at this 
early stage. The preliminary tests presented here aim at 
demonstrating the ability of the method to yield signifi- 
cant constraints on a one parameter family for the form 
of 5* as given below. We test the method using catalogs 
of mock data of dynamically relaxed tracers. The system 
is assumed to exist in a galaxy made of a spherical dark 
matter halo and a baryonic disk, neglecting the gravita- 
tional effects of a bulge . The grav itational potenti al of 
the h alo is taken as (|Sakamoto et al.,,2003; .Dinescu et al.l 
I1999D 



v,Hn [l + {r/df] 

2 3 

r r ,2+^2 ; 



for r < r 
otherwise , 



(3) 

where ^"0 is a constant ensuring the continuity of the 
force field per unit mass, —'^'^t.^i^, a,t r — r^^^. We take 
y^, = 128km s~\ d = 12 kpc and r^^^ — 170 kpc as in 
ISakamoto et al.l ( 2003| ). For the potential of the disk we 
work with the form (Miyamoto & Nagai 1975), 



^ disk 



-GAL 



(4) 



where d = 12 kpc, M,^,.,. = 10"M (n, a = 6.5 kpc and 
b = 0.26 kpc (|Dinescu et al.lll999f ). The mock observa- 
tions are generated as follows. Particles are placed in 
the halo by a random Poisson sampling of an underling 
number density n{r) cx 1/r^. The corresponding veloc- 
ities are selected randomly from a gaussian distribution 
with zero mean and a r.m.s proportional to the local es- 
cape speed such that the total kinetic energy is half the 
absolute value of the potential energy, as implied by the 
virial theorem. The equations of motion of the system 
are then solved numerically for 10 Gyrs to ensure dynam- 
ical relaxation in the gravitational field given in Eqs.[3]& 
m Outputs of particle positions and velocities obtained 
from this procedure are then identified as the mock cat- 
alog to be used for testing the method. As our model 
gravitational potential we use the forms in Eqs. [3] & S] 
but with M^j^j^ and v^^ ^ scaled by a factor a according to 
{M^.^^,v^^^) (aM^.^j^, aWj^^). The tests are performed 
for one free parameter, a. Following the scheme above, 
the mock catalog provide the initial conditions at time 
tg for moving the particles for another 10 Gyrs for an 
assumed value of a. Snapshots of particle positions are 
then tabulated at 1000 uniformly distributed time steps 



between and -I- 10 Gyrs. This is done for 17 val- 
ues of a spanning the range 0.7 — 1.5 linearly. The tests 
will demonstrate that the method is able to recover the 
correct value, a = 1, within an acceptable uncertainty. 

3.1. Tests vi^ith zero measurement errors 

We start with the statistic N computed from snap- 
shots corresponding to an assumed a. The filled circles 
in the left panel of Fig. [T] show TV as a function of a, for 
three values of the number of tracers, n = 200, 300 and 
500. The error-bars represent the r.m.s scatter, a. The 
horizontal solid lines indicate the limit, N = n/2, corre- 
sponding to the three values of n. For a close to unity, 
N is indeed near n/2 even for n = 200. For better quan- 
tification of these results we plot in the right panel the 
quantity x^(a) = [N — n/2)^/(T^. Curves of become 
narrower as n is increased, but already with n — 200 
significant constraints on a could be derived. Note that 
< (TVs - n/2)2 > /cr^ = 1 -I- so that x^ = 3 correspond 
to a 2tT deviation of Ns- We see that n — 200 already 
constrains a to better than 25% at the 2a level. 

We now turn to the distribution function of r and ^ 
as described in 52 We define the cumulative distribu- 
tion function (hereafter CPDF), P(< j/i), at any time, 
as the fraction of particles having y < j/i, where y is 
either r or ^. Hereafter, we will treat the CPDF ob- 
tained from the ensemble of 1000 snapshots as the un- 
derlying actual model CPDF and denote it by P^. In 
the left panel of Fig. [2 we illustrate the differences be- 
tween the various CPDFs obtained with 300 tracers. The 
thick solid curve is the CPDF P(< r) obtained from 
the "observed" positions in the mock catalog, while the 
dashed, thin solid, and dot-dashed lines correspond to 
P„(< r) for a = 0.7, 1, and 1.5, respectively. The CPDF 
P^(< r) is sensitive to a. The value a = 0.7 broadens 
the distribution of particles towards larger radii, relative 
to the "observed" CPDF, while a = 1.5 concentrates the 
particles nearer to the center. The CPDF P„(< r) for 
a = 1 (smooth thin solid line) seems consistent with the 
"observed" CPDF (thick sohd). To quantify the differ- 
ences between model and observed "observed" CPDFs, 
we compute the Kolmogorov-Smirnov (KS) "distance" 
Da = sup |-P-„(< r) — P{< r)\ where the dependence on 
a is only through P^^{< r). Given Da and the number 
of tracers, n, we compute the significance level, Q^s' 
the hypothesis that P{< r) and P„(< r) represent the 
same distribution. In the right panel of Fig. [2 we plot 
Qj^g as a function of a for n = 200, 300, 500 and 800, as 
indicated in the figure. All curves peak near the value, 
a = 1. The only uncertainty in determining a is due 
to the finite number of tracers, hence the peaks are nar- 
rower for larger n. For all values of n considered here. 
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the choice a = 1 is never "rejected" at more than the 
32% confidence level (i.e. Q^^ < 0.32), equivalent to a 
ItJ level for a normal distribution. Thus the method pro- 
duces unbiased estimates of a within the la errors. For 
n = 200, deviations of a larger than ^ 30% from unity 
are rejected at more than the 2a level (i.e. Qj^g < 0.045). 

We now explore the CPDFs of ^. The CPDF P(< C) 
computed from the "observed" particle positions in the 
mock catalog depend on the assumed a through the peri- 
centers, r^^i^, and apocenters, r^,^^, of the orbits. This is 
in contrast to the observed P{< r) which is completely 
independent of a. The three thick lines in the left panel 
of Fig. [3] are the CPDFs P(< ^), computed from the 
"observed" positions of 300 particles, for a — 0.7, 1 
and 1.5, as indicated in the figure. The three nearly 
overlapping thin dashed, thin solid and thin dot-dashed 
lines show P„(< ^), respectively, for these three a val- 
ues. The dependence ol P^{< £,) on a is very weak, in 
contrast to P^{< r) which is very sensitive to the as- 
sumed a. This weak dependence is not entirely unex- 
pected since ^ eliminate the overall length scale of the 
problem. For a = 0.7, the difference between P„(< £,) 
(thin dashed line) and P{< £,) (thick dashed) implies that 
most of the "observed" positions are closer to the peri- 
centers than the positions in the snapshots. For a = 1.5 
(thin and thick dot-dashed lines), the "observed" posi- 
tions are closer to the apocenters than the positions in 
the snapshots. For the correct value a = 1, the distri- 
butions P^{< C) (thin solid) and P(< ^) (thick solid) 
appear to be consistent. To test whether the difference 
between P^{< £,) and P{< £,) is large enough to rule an 
assumed a we show in Fig. [3] the corresponding quantity 
Qj^g as a function of a for n ~ 200, 300, 500 and 800, as 
indicated in the figure. As is the case for the CPDFs of 
r, the value a = 1 is never "rejected" at more than the 
la level. All curves are narrower than the correspond- 
ing curves in the right panel of Fig. [21 implying a better 
ability to constrain a with the distribution of ^ than of 
r. For n = 200, deviations of more than 20% from a = 1 
are rejected at more than the 2a (Q^s < 0.045) level. 
For a given confidence level, the range of a (i.e. con- 
fidence interval) constrained by the distribution of ^ is 
narrower by a factor of 2 than the range constrained by 
distribution of r. 

The statistic \W/2T\ computed form the mock obser- 
vations with n — 200 is presented as a function of a by 
the filled circles in the left panel of Fig. [D Here T is 
computed with the "observed" velocities and the depen- 
dence on a is entirely due to the W term (see Eq. 
For each a, the error-bars show the r.m.s scatter, a, of 
\W/2T\ estimated from snapshots generated with that 
a. As a consistency check we confirm that the departure 



from unity of the ratio W/2T, obtained by averaging 
W and T over the snapshots, is completely insignificant 
compared to the la scatter. The deviation from unity of 
\W/2T\ computed from the mock observations is shown 
in the right panel of the Fig.Has = [\W/2T\ - l]^/a^. 
For n = 200, this statistic constrains a to better than 
15% at the 2a confidence level. This is better than the 
N statistic which constrains a to better than 25% at the 
same confidence level. It also fares slightly better than 
the ^ distribution. 

3.2. Effect of measurement errors 

So far we have assumed zero errors in the phase space 
coordinates. We present here only partial tests of the 
method when applied on noisy data. The amplitude of 
the errors depends on the distance of tracers from an 
observer at the solar position rather than the galacto- 
centric distances. We assign a 15% error in the distance 
of a tracer from an observer present in the Galactic disk 
at 8kpc from the Galactic center. For GAIA a 10% par- 
allax distance error corresponds to objects of 15 mag at 
lOkpc from the observer. Errors in parallax distances 
scale quadratically with true distance and become very 
large at distances of tens of kpcs even for GAIA and SIM. 
Here we assume that distances of far away tracers are de- 
termined by other means so that a linear scaling of the 
errors is maintained. We also perturb the radial velocities 
and proper motions with errors that scale linearly with 
distance. The amplitude velocity errors is assumed equal 
in all three directions and is normalized to a r.m.s value of 
lOkms"^ at a distance of 20 kpc. For comparison, at 15 
mag, GAIA after 5 years of operation will provide radial 
motions within an accuracy of 10— 15kms~^ and proper 
motions within an accuracy of 1km s^^ at a 20 kpc. 

We do not show results for P(< r) computed with er- 
rors included as it yields less significant constraints than 
the remaining statistics. The results for N and ^ are 
shown in Fig. [5] and Fig. [6l respectively. For n = 200 
and 300, there is little difference between the curves of 
and Q^g in these figures and in Figs. [1] & [3] corre- 
sponding to zero errors. For larger n the effect of the 
errors is more pronounced, especially in Qj^g . The reason 
for this behavior is that sampling errors resulting from 
the finite number of tracers are dominant over measure- 
ment errors for the smaller values of n. As n increases, 
measurement errors become more pronounced resulting 
in the significant reduction of the value of Q^g at a = 1 
for n = 800. The N statistic seems to be more resilient 
to measurement errors than the ^ distribution. The virial 
ratio \W/2T\ computed with noisy mock data is shown 
in Fig. [7] with the same notation as Fig. d] correspond- 
ing to results with zero measurement errors. In comput- 
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ing T, the mean of non-vanishing quadratic terms due 
to velocity errors have been removed. Measurement er- 
rors seem to affect the ratio more than the other statis- 
tics. Still, this ratio constraints a slightly better than 
the other statistics. 

Overall, measurement errors of the amplitude we con- 
sider here have not degraded the ability of the method at 
constraining a. Therefore, the method does not require 
unrealistically accurate data. 

3.3. Effect of gradually growing disk 

The method assumes a constant gravitational poten- 
tial. However, the Galactic disk may have grown sub- 
stantially in the last 8Gyr or so. Here we check whether 
this gradual growth seriously hampers the application of 
the method. We solve for the orbits of tracers assum- 
ing that the mass distribution of the disk grows like t'^^^. 
The method is then applied to the resultant distribution 
of tracers in phase space assuming a constant disk. In 
Fig. [8] we show the confidence level Q^^ as a function of 
a where in this case the mass in the disk is assumed to 
be known and a describes the ratio of the assumed to 
the actual value. 

4. Concluding Remarks 

The ergodic method presented here requires a para- 
metric functional form for the Gravitational potential, 
but does not impose any special symmetry on the mass 
distribution. The method assumes that the observations 
of a class of tracers are spatially complete. Tracers with 
observed distances smaller than do, could be present be- 
yond do in snapshots at later times. Therefore, Ob- 
servational selection against tracers at distances > do 
will make the statistical comparison between observa- 
tions and snapshots at later times extremely difficult. 
The completeness is not too demanding a condition for 
tracers like globular clusters and Galactic satellites for 
which future observations should be accurate enough for 
quite large Galactic distances. 

We have presented only partial testing of the method, 
with only a one parameter family for the form of the 
Galactic potential. Method is able to constrain the pa- 
rameter to a good accuracy with measurements errors 
that are even larger than those expected to be achieved 
by future data. The tests show that the method could 
provide unbiased constraints on the Gravitational poten- 
tial, but a more elaborate testing which includes a more 
realistic treatment of the errors should be done. Observa- 
tions will likely assign distance and velocity measurement 
errors to tracers on an individual basis. Therefore, ran- 
dom errors and systematic biases tailored to the specific 



sample of tracers used by the method could be deter- 
mined robustly. 

The accuracy of the method is mainly limited by the 
number of tracers. The most obvious tracers are globular 
clusters and Galactic satellites. Our Galaxy includes 158 
known globular clust ers, and 23 known satellites (e.g. 
ISimon fc Geh^ l2007^ . Distance measurements of RR 
Lyrae stars from their period-luminosity relation will be 
greatly improved by GAIA and SIM calibration of the 
zero point using a nearby sample of these stars. There- 
fore, luminous halo RR Lyrae stars could significantly 
enlarge the sample of tracers. 

The method requires a system of tracers in dynam- 
ical equilibrium in the current Galactic potential. A 
pre-requist for dynamical equilibrium is that any recent 
changes in the Galactic potential must have occurred on 
a time scale longer that the dynamical time of the system 
of tracers. Spectroscopic studies stars in the Galaxy do 
not pr esent evidence for subst antial mergers in the last 
SGyr (|Gilmore et all [20021 : iHelmi et all [20061 ). implying 
a nearly static gravitational potential. We have demon- 
strated that adiabatic growth of the Galactic disk is not 
expected to pose a problem for the implementation of the 
method. The effect of any deviation from that should be 
modeled. 

One issue which is exclusive to using disk stars as trac- 
ers is the effect of tra nsient perturbations on the dynam - 
ics of those stars (e.g. iFamaev et al. |[200l[Quiiie3[200i). 
However, the transient effects cause velocity perturba- 
tions at the level of 10 km/s which is much smaller that 
the total velocities of tracers. Observational uncertain- 
ties are larger than that and do not seem to cause sig- 
nificant biases in the results as indicated by the tests 
described above. Another issue to be considered is halo 
substructure which, in principle, could act as a stochastic 
component in the gravitational potential. Such a compo- 
nent is extremely difhcult to model in the method pro- 
posed here. We offer the following argument demonstrat- 
ing that substructure should not have an important effect 
on the long term dynamics of tracers. In the limit of fast 
encounters, a tracer passing a substructure at distance 
b will change its velocity by Vi « g{b)b/V where V is 
the relative velocity and g « GM/b^ is the gravitational 
force field of the substructure assuming that b is larger 
than its tidal radius. Performing the usual summing in 
quadratures over encounters occurring in one orbital time 
we get that the r.m.s change < Vi >^/^« MV^R/V 
where n is the number density of substructures and R is 
the distance travelled by the tracer in one orbital time. 
Taking = GM^,„{R)/R with M^,„{R) the mass of 
the Galactic halo within radius R, we get the condi- 
tion ~ M^^VnR? > M„„ for having < >i/2^ V. 
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This means that neither single encounters nor collec- 
tive stochastic effects can dominate the long term evo- 
lution of tracers even if the fraction of mass in substruc- 
tures is large w hich is contrary to recent simulations (e.g. 
IColombill2008l ) which show that the smooth component 
greatly dominates the mass of Galactic size halos. 

When details of future data from SIM and GAIA be- 
come available, all robust information about the distri- 
bution of b aryons in the Galaxy should be used (e.g. 
Ir^ in et al.li2003 ) in order to place tight constraints on 
the Galactic dark matter. The validity of the method 
should be tested with mock data that match the obser- 
vations as much as possible and with the best possible 
available Galactic models. 
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Fig. 1. — Left: The statistic N as a, function of a shown as the filled circles with error-bars representing the r.m.s scatter, 
a (see text). Three values of the number of tracers are considered, n — 200, 300 and 500 and the three horizontal lines 
indicate corresponding limiting values = n/2. Right: Curves of = (iV — n/2)^/a^ as a function of a for several 
values of n, as indicated in the figure. 




Fig. 2. Left: CPDFs of the distances for 300 tracers. The thick line shows P{< r) computed from the observed 
"distances: at t = t^, while the three smooth thin lines represent the model CPDF P^{< r) for three values of a, as 
indicated in the figure. Right: The confidence level Qks that P(< r) and P^{< r) represent the same distribution, as a 
function of a. 
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Fig. 3. — Left: CPDFs of ^ for 300 tracers. The three thick hnes represent "observed" P(< ^) computed with r^.^ and 
r^^^ corresponding to a = 0.7, 1 and 1.5, as indicated in the figure. The nearly overlapping three smooth thin lines 
represent the model P^{< C) corresponding to the same choice of a values. Right: The confidence level Q^^ as a function 
of a. 




Fig. 4. — Left: The filled circles show the ratio, |VK/2r|, versus a from the observed distances and velocities of 200 mock 
tracers at tg. The error-bars show the r.m.s. scatter, a, computed from snapshots a.t t > t^. The horizontal line indicates 
the virial theorem value for W/2T. Right: The quantity = [|W^/2T| - 1]^/(T^ versus a. 



8 



300 




Fig. 5. — The same as Fig. [l] but with degraded mock observations by random errors as described in the text. 




Fig. 6. — The same as the right panel of Fig. [3l but with degraded mock observations. 
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Fig. 7. — The same as Fig. 21 but with degraded mock observations. 
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Fig. 8. — The confidence level Q^g as a function of a, 
assuming a gradually growing disk (see text). 
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